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Abstract 



We study effect of noise on the tunneling process of the magnetization 
in a sweeping field from the viewpoint of the nonadiabatic transition. The 
tunneling probability depends on the properties of the noise, e.g. the direc- 
tion, the amplitude, and the relaxation time of the autocorrelation function. 
We investigate the tunneling process in the presence of noise by solving the 
time-dependent Schrodinger equation numerically, and analyze its parame- 
ter dependence. We estimate the effective tunneling gap making use of the 
Landau-Zener-Stiickelberg (LZS) formula. The dependence of the gap on the 
sweeping velocity reproduces the corresponding experimental data qualita- 
tively. For fast sweeping field the estimated gap has no velocity dependence 
and this gap seems to give the true gap. We, however, show that a correc- 
tion for the estimated energy gap is necessary and we propose an correction 
formula. 
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I. INTRODUCTION 



The quantum tunneling of the magnetization (QTM) has become one of the new topics 
in the condensed matter physics. At low enough temperatures phenomena originating quan- 
tum dynamics can be observed in mesoscopic and microscopic systems. Recently quantum 
dynamical effects have really been observed in the hysteresis phenomena of the nanoscale 
molecular magnets such as Mn 12 , Fe 8 , and V15, and the analyses of these phenomena are in 
progress [|T| [3~9f| • These mesoscopic quantum phenomena have been actively investigated not 
only for academic interests but also for future applications, e.g., in quantum computers. 

We have studied this tunneling dynamics from the viewpoint of the nonadiabatic tran- 



sition |20|-Pq1. The transition probability for the nonadiabatic transition is given by 



the Landau-Zener-Stiickelberg (LZS) formula in the pure quantum systems, i.e., without 
noise [27|. We have shown that the LZS formula is useful for the estimation of the energy 
gap at the level crossing point. 

In the realistic system the effect of the coupling to the environment plays a crucial 
role in such phenomena. At low enough temperatures, below thermally activated regime, 
phonon-mediated relaxation can be neglected and the effects of the ubiquitous hyperfine 
field within the molecule and the dipole field of the other molecules become significant 
for the relaxation process. They are considered to be the main source of decoherence. 
The 'nonexponential relaxation' phenomena is one of the typical examples. Time decay 
of the remanent magnetization measured at low temperatures at the early stage has been 
reported not to obey usual exponential function but shows a square-root function both in 
the experimental and in the theoretical analyses p9|-|3T||. 

The LZS tunneling probability is a function of the sweeping velocity and the energy gap 
at the resonant point. Because the sweeping velocity and the transition probability can 
be obtained in the experiments, the gap can be estimated from the LZS formula. In the 
experiments of Wernsdorfer and co-workers, the tunnel gap has been estimated for molecular 



magnet Fes The estimated energy gap shows a dependence on the sweeping rate 

due to the noise effects, which does not take place in the pure system. The estimated gap, 
however, is almost independent of the sweeping velocity in the range of fast sweeping velocity 
while it decreases when the velocity decreases. The estimated gap in the fast sweeping range 
seems to give a true tunneling energy gap, judging from the fact that the velocity dependence 
does not exist in the pure quantum system. Namely the noise effect seems to be negligible in 
the fast sweeping range. We clarify that this is not always true and show that a correction 
of the estimated gap is necessary. We propose the relation between the true value of the 
gap and the estimated one. 

In this paper we provide a stochastic noise field and investigate the time evolution of the 
system by solving the Schrodinger equation numerically. Kayanuma has given the analytical 
formulas of the transition probability in the stochastic noise for special cases, and also given 



numerical results obtained by Wiener-Hermite expansion [32|. The effect of the boson field 



has been also studied and the transition probabilities have been given analytically for various 
parameter range [Q. Thus far several studies have been performed on the noise effect where 
the noise is parallel to external field. As far as we know, however, the effect of noise applied 
in other directions has been rarely dealt with and also been limited to the special case. Here 
we also examine more general cases and investigate direction dependence of the noise. 
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From the averaged time evolution we study how the transition probability, as a function 
of the sweeping velocity, depends on the direction, the relaxation time, and the amplitude 
of noise. We also compare our results with the analytical results for the limiting case |32| . 
Furthermore, we investigate the relationship between the tunneling gap estimated from fast 
sweeping rate and the true value, and propose a correction form of the estimated energy. 

The layout of this paper is as follows. In Sec. ||, we explain the method used in this 
paper. In Sec. |TT|, we study the dependence of the transition probability on the damping 
factor and the amplitude of noise. In Sec. |TV| , we investigate sweeping velocity dependence 
of the energy gap estimated from the transition probability and propose the correction of 
the estimated gap by LZS formula. A summary and discussions are given in Sec. [V] 



II. METHOD 

Here we consider a two level system in a sweeping field with noise. The Hamiltonian is 
given by 

H{t) = 2TS X - H(t)S z + H noisc {t), (1) 

where H(t) is the sweeping field H(t) = (—H + ct) and TC noise (t) describes the noise. We 
investigate the transition probability under the noise applied in the x-, y-, or z-direction, 

H noise {t) = h x {t)S x , h y (t)S y or h z (t)S z . (2) 

We solve the time evolution of the time-dependent Schrodinger equation (TDSE) 

ih^Mt)) = n(tmt)h (3) 

where denotes the wave function of the spin system at time t. Hereafter we set 

h = 1. In this time evolution, we make the exponential decomposition exp(— iTi/St) ~ 
exp(— iH x At) exp(— iTC y At) exp(— iH z At), where TC X depends only on the x-direction op- 
erator and similarly in the y- or z-direction. We use the rotating operator U(a — > b) 
(a, b = x,y, z) which rotates the quantum axis from a-axis to 6-axis and make each exponen- 
tial diagonalized. Using this method ||34|| , we can also simulate the case that the noise has 
components in all directions (h x , h y , h z ). The present system is a two-level system whose 
adiabatic energy levels are shown in Fig. [I] (7i noisc (t) = 0) as a function of the magnetic field 
H. There the avoided level crossing occurs at H = 0. 

As the initial state, we set H(t = 0) = — H < 0, and put the system in the ground state 
for this field. For large H > 0, the initial state is almost the down state (M z = —1/2 ), 
i.e, \ijj(t = 0)) = |— ). Without noise the probability to remain in the ground state after the 
crossing is given by the Landau-Zener-Stiickelberg (LZS) |27j formula, 

p = \(+Mt = °°)}\ 2 (4) 
= 1 - exp J , (5) 

which is a function of the sweeping velocity and the energy gap (AE = 2T) at the avoided 
level crossing point. 
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When noise is present, it disturbs the quantum process and hence this formula for the 
probability changes. We will investigate how the probability depends on the noise. We 
denote the transition probability under the noise field in the x-, y- and ^-direction by p x , p y 
and p z , respectively. 

It is naturally assumed that noise sources such as the hyperfine field and the dipole 
filed independently cause a change of the internal field with some restoring force on the 
target spin. The time evolution of this random noise can be given by a Langevin equation 
(Ornstein-Ulenbeck process), 

/ i (t) = - T / i (t)+r 7 (t), (6) 

where r)(t) is a white gaussian noise 

( V (t))=0 and ( V (0) V (t)) = A 2 5(t), (7) 

and 7 is the damping factor due to the restoring force. 

This random noise h(t) has a gaussian distribution and an exponential-decaying auto- 
correlation function. 

(h(t)) = and (h(O)h(t)) = ^exp(--), (8) 



where r = I/7 is the relaxation time. Here \jA 2 /2^ is regarded as the amplitude of the 

noise and abbreviated by B = J ' A 2 ~/2y. We prepare h x (t), h y (t), and h z (t) independently. 

In realistic situations, the noise has components in all directions. We also study this 
case, where we let the noise have the same amplitude as in the case of one direction, that is, 
(h 2 + hy + hi) = B 2 . We take an average over 1000-4000 samples to estimate the physical 
properties. 



III. DEPENDENCE OF THE TRANSITION PROBABILITY ON THE NOISE 

In this section we investigate the dependence of the transition probability on the damping 
factor and also on the amplitude of the noise. The dependence of the noise direction is also 
investigated. 

Here we focus on the case of slow sweeping velocity because the dependence on the 
direction, the amplitude, and also 7 appears more clearly than in the fast sweeping case 
(see next section). The field is swept from H(t = 0)=— 5.0 to H(t = £fi na i)=5.0 for the 
system of T = 0.02. Because H(t = 0) ^> T, the initial state is down spin state. We set 
the sweeping velocity to be very slow (c=0.0005), such that the system almost shows the 
adiabatic transition, i.e., p ~ 1 when B = 0. 

Figures |2] (a) and (b) show the time dependence of the noise of 5=0.02 at the values of 
7 = 0.001 and 7 = 0.32, respectively. These represent the typical slow and rapid fluctuations. 
Figure |3| (a) shows dependence of the transition probability on the damping factor 7 keeping 
the amplitude constant, 5=0.02, and Fig. |3] (b) shows that for the case of 5=0.005. In 
the former parameter set, the energy gap without noise (2r) and the noise amplitude (B) 
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have the same order of magnitude while in the latter parameter set, the noise amplitude B 
is much smaller than the gap. 

As we naturally expect, the noise effect is larger for larger amplitude B of the noise. 
Although the absolute values of the transition probabilities are different, they have a similar 
dependence on the damping factor. The noise of x- or ^/-direction has stronger influence 
than the noise in z-direction. This is because the transverse field causes the transition to 
the other level directly. For very small value of 7, the noise can be considered to be static 
and p x , p y , and p z are approaching to a constant as 7 — > 0. The probabilities p x , p y , and 
p z decrease rapidly around 7 = 0.02. This value corresponds to the inverse of the transition 
time, 

_ = _ = 0.025. (9) 

ttv - 

The transition time is the time scale to pass the avoided cross region. This indicates that 
when the time scale of noise fluctuation becomes shorter than that of the transition time, 
the transition probability changes very sensitively regardless of the noise amplitude. The 
probabilities p x and p y take the minimum values at 7 = 0.1 — 1.0 while p z has the minimum 
value at 7 = 0.02 — 0.05. The minimum values of p x and p y are close to each other and they 
are smaller than that of p z . When 7 becomes very large, the effect of noise becomes small 
due to a kind of motional narrowing. Let us study these dependences in detail. 



A. Slow fluctuation limit 7 — > 

In the region of small 7, i.e. in the slow fluctuation case, there is a distinguishable 
difference between p x and p y , while p y and p z approach to each other. This tendency is 
found in both figures ^| (a) and (b). These dependences can be understood by an analysis 
of the static distribution of noise. 

When the noise is applied in the x-direction the magnitude of transverse field is T + h/2. 
Thus, p x is obtained by averaging over the distribution of h: 

Px = (plzs (r - r + h/2)) = (1 - ex P (- 27r(r + V2)2 ) ^ , (10) 

where (• • •) denotes 

1 r°° ( h 2 \ 

Here the notation (T — > X) means the replacement of V by X in the LZS transition probabil- 
ity Plzs- In the same way, when the noise is applied in the ^-direction, the total magnitude 
of transverse field is \Jt 2 + (h/2) 2 , and we have 

^( te (r.^)).(i- esp (-M))). (12) 

When the noise is applied in z-direction, the field shifts from H(t) to H(t) + h and has no 
effect on the transition probability, 
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V* = (Rzs) = U - exp ] \ . (13) 



c 

Explicitly evaluating the average with Eq . (^ [T0| , [T2| , p~3D , we have 



1 / 1 /7T„\ 2 27T^ 2 



and 



^ = 1 -V2^ eXP (-^-)' (15) 



/ 2vrr 2 \ , , 

p z = 1 - exp , (16) 



where we define 

When T = 0, the transition probability p x (= p y ) has been obtained analytically in the 



limit of 7 — > calculating by the infinite series expansion |32fl. This p x coincides exactly 
with the special case (T = 0) of (Eq. |I4]) , 

fc -(i--p(-^)>-i-(i + ^)-*. 

Here we find that p y > p z and p y > p^, because 

Ut>D = 2&(± + (19) 

and 



1 / 2ttT 2 \ f ( 1 



71", 



2^ 



^ = V5j#u exp {-—) r p u j - i ) > °- (20) 

For the case (a) (r = 0.02, B = 0.02, c = 0.0005), we have p x = 0.8724, p y = 0.9965, and 
p z = 0.9934, and for the case (b) (r = 0.02, B = 0.005, c = 0.0005), we have p x = 0.9879, 
p y = 0.9939, and p z = 0.9934. As shown in Figs. |3| (a) and (b) the values obtained by 
simulations approach these values. 



B. Fast fluctuation limit 7 — ► 00 

As 7 becomes very large, the noise in Eq. (§) becomes a kind of the white noise 

(h(0)h(t)) = —6{t) (7 - 00), (21) 
7 
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but its amplitude becomes zero 



2 D2 

lim / 5{t)dt = ^ l B 2 -> 0. (22) 



7— >oo 
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Therefore the effect of noise becomes small and the probabilities increase again and reach the 
probability of LZS regardless of the direction of noise. This effect resembles to the motional 
narrowing which often appears in experiments of spin resonance. Namely, the change by the 
noise is averaged out. 

On the other hand, under the condition that B 2 ^> 1, p z in the limit 7 — > 00 has been 
given analytically by Kayanuma 

pi= l-exp(-4,rVc) {a) 
This formula yields p z = 0.500 for the parameters T = 0.02 and c = 0.0005. Thus the value 



of p z is expected to change from the value of LZS to 1/2 when B = J A 2 /2 / y increases. Let 
us study how the probability changes as A becomes large for large 7(=100). A-dependence 
of p z is shown in Fig. 0], where p z approach 0.5 as A becomes large as we expect. The fact 
that p z = 0.5 suggest that when the noise is large, large energy pours into the system during 
the crossing and the distribution over the ground state and the excited state equalizes. 
These situations are similar to those with thermal bath at very high temperatures. We also 
investigate A-dependence of p x and p y (Fig. 01) and we find p x and p y approach 0.5 more 
rapidly than p z as A becomes large. 

For the special parameter (r = 0), p x in the limit 7 — > 00 has been also given analytically 
by Kayanuma fl32fl , 

1 — exp(— -kB 2 /c) , , 

P*(r = 0) = ^ '-J-. (24) 

If we perform a simulation for large values of 7, p x does not follows this formula, but it 
approaches 0. For example we find p x (T = 0) ~ for (B = 0.02, c = 0.0005, and 7 = 100), 
while the above formula gives p x (T = 0) = 0.459. The numerical result is due to the 
abovementioned narrowing effect. Where does this inconsistency come from? In analytical 
treatment, the field is swept from —00 to 00, while in the simulation the field is swept in 
a finite region, i.e., from — H to H . If we set range of sweeping field to be finite in the 
Kayanuma's derivation, it is found that p x (T = 0) = [55|. Thus the limits of 7 —>■ 00 



and H — > 00 are not exchangeable. A detailed analysis of these limiting cases will be given 
elsewhere. 



IV. SWEEPING VELOCITY DEPENDENCE OF THE EFFECTIVE GAP 

As we mentioned in the introduction, the sweeping velocity and the transition probability 
are observed experimentally and the gap is estimated by making use of the LZS formula. 
In this section, we investigate the dependence of the estimated energy gap (AE) on the 
sweeping velocity. 
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In the absence of noise, the gap at the crossing is given by AEq = 2T. Similarly in the 
case with noise, the effective gap can be estimated using LZS formula as 




AE eS = \l--\n(l-p). (25) 
Let us define the following ratio, 



which indicates the strength of the noise effect. If this ratio is nearly equal to 1, the effect of 
noise is negligible. Here we investigate how the noise effect is changed due to the sweeping 
velocity. 

We show the sweeping rate dependence of the effective gap AE efi for various damping 
factors in Fig. [S]. Here we set the parameters B = 0.02. Figures (a)-(c) show the cases 
when the noise is applied in x-, y-, and z-direction, respectively. Figure || (d) shows the 
case when the noise has components in all directions (h x , h y , h z ). Here we let the noise have 
the same amplitude as in the one-direction case, i.e. (h x + h 2 + hi) = B 2 . The following 
tendencies appear for the cases (a)-(d). For large sweeping velocity, the effective gap is found 
to be constant. In the range of slow sweeping velocity, the effective gap deviates from this 
constant and decreases roughly linearly with the logarithm of the sweeping velocity. These 



results are consistent with the experimental results []n|fuj, where the estimated gap AE e g is 
constant at larger velocity and exhibits a deviation in the low sweeping velocity range. We 
find that this deviation also significantly depends on the damping factor. As the sweeping 
velocity gets slow, the system stays near the avoided crossing point for a long time, where 
the adiabatic motion is largely disturbed, and then effectively the gap becomes small. 

Let us consider the dependence of AE c s / AE on c and 7. The field crosses the avoided 
level crossing region in a time of the order of r tr = Y/c. When 7 is much smaller than 
r t 7 = 7tr, the noise can be regarded to be static. We can estimate AE e g/AE in the 
static limit by using the static noise results for p x , p y , and p z in Eqs. fllil) , (|15|) , and (|I6|). 



Substituting them into Eq. (^), we have the velocity dependence of AE^ S / AE (d = x, y, 
or z). These dependences are shown in Fig. ^| and also shown by the solid lines in Fig. |5| 
(a)-(c). We find that these solid lines agree well with the data for large sweeping velocity c. 
The shape of the solid line in each direction is different from one another. 

When 7 increases, the deviation of the effective gap from the solid line becomes large in 
the slow sweeping region and the critical value of the velocity, below which the AE^ S / AE 
(d = x, y, or z) deviates, increases. The effect of noise is found to be less important for 
z-direction. The data of 7 = 0.32 and 1.0 in Fig. |5| (c) are exceptional and these points 
represent the behavior in the narrowing region in Fig. [3], that is, in the region where the 
transition probability recovers. Thus the noise is inefficient there and A e g behaves differently. 

At large c, r tr = 7^ is small thus j tT is large. Therefore, even the noise with large 7s 
can be regarded to be static and these data still follow the solid line. If we assume that the 
noise in the realistic system has a certain finite relaxation time, the transition probability 
corresponds to that of static limit case for very fast sweeping. 

We know that AE c g has no sweeping velocity dependence for fast sweeping but its value 
can be different from the true value of the gap. Let us estimate the correction of the 



S 



gap. When the sweeping velocity is large enough, i.e. 2irT 2 jc <C 1, the LZS probability is 
approximately given by 

2nT 2 7r o 

PLZS - — = Y C AE °- (27) 

Taking the same limit (7 — > 0) in the formulas (0), (0), and (|I6|), we have 

» = (^(r^r + £)) = ?^^ = £(A^, (28 ) 



ii 



P, = (plzs (r - A 2 + (V2) 2 ) ) = ^ 4 J = ^(A^ff) 2 , (29) 



and 



v 2vrr 2 , . 

Pz = <Plzs> = • (30) 



Thus from p z we can obtain the true value of the gap AE = 2T. On the other hand, 
from Eqs. (|27|) and (^8)), we have 



AE^\ 2 . B 2 



AE J 1 + 4P- (31) 

Probability p y has the same correction. For the present case B 2 /T 2 = 1, we have 
AE^ S /AE =AE^ S /AE =1.118, and AE* S /AE =1, which agree with the simulation results 
(Fig. l(a)-(c)). 

For the isotropic static noise, transition probability is given by 

_ p x +Py+Pz 27rr 2 B 2 

P = 3 = — + Tc- (32) 



Then we have 



p fAE cS \ 2 , B 2 

1 + 7^, (33) 



that is 



AE e{i = ^l + -^AE . (34) 

Thus we find that the true gap AEq is obtained within a correction of 0(j&). Even if 
the noise amplitude is large, the estimated value is proportional to the true value as far 
as the properties of noise are not changed. In the present case, because B 2 /T 2 = 1, 
AEgg/ AE — 1.083, which also agree with the simulation results Fig. |5| (d). The effec- 
tive gaps are always larger than the true value and these give the upper limit of the true 
gap. 



9 



V. SUMMARY AND DISCUSSION 



We have studied the effect of noise on the nonadiabatic transition in the sweeping field. 
The transition probability depends on the amplitude and the relaxation time of the noise. 
The probability also depends on the direction of the noise. In the slow fluctuation limit of 
the noise, we can regard the noise to be static, where we obtain the analytical formulas of 
p x , p y , and p z , and find the relations p y > p z = p^zs an d p y > p x - When the time scale of 
the noise fluctuation comes close to that of the transition time, the transition probability 
drastically changes. When 7 increases further, the probabilities approaches that of LZS. 
There the time scale of the noise fluctuation is much smaller than the transition time. We 
found that there exists an important correlation between these two time scales and the effect 
of noise. 

In the regime of slow sweeping, the estimated gap as a function of the sweeping velocity 
depends on the direction and the damping factor of the noise. When the sweeping field 
becomes large, the estimated gap gets constant regardless of the damping factor of the 
noise. It seems that the noise has little effect in this regime. However the effect of noise 
does not disappear. These constant values depend on the direction of the noise, AE^ S = 
AE^ S > A£g ff = AE true . When the amplitude of noise is not small, the noise effect is not 
negligible. We obtained the relationship between the estimated gap and the true one. 

In the Wernsdorfer and co-worker's experiments, Fe§ clusters show a little difference of 
the estimated gaps between the isotropic modified samples [E^jnj. These differences may 
come from the abovementioned noise effect. In order to clarify the situation, it is necessary 
to deal with the realistic model (S = 10) and to estimate the gap for M z = ±10, where the 
effect of noise in the x- and y-direction would be a little different from the present case. 

We believe that our study captures an essence of the noise effect in the relaxation process 
of spin dynamics and gives important knowledge to following studies. We will investigate 
details of a more realistic model in the near future. We hope that this work provides a basic 
information for the manipulation of microscopic or nanoscale magnetic devices in the future. 
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FIGURES 



FIG. 1. Avoided level crossing. 

FIG. 2. Time dependence of the noise, (a) the slow fluctuation case (7 = 0.001) and (b) the 
rapid fluctuation case (7 = 0.32). 

FIG. 3. 7 dependence of the transition probability in the case of slow sweeping velocity, (a) 
5=0.02 (b) 5=0.005 

FIG. 4. yl-dependence of p x , p y , and p z for large 7(= 100). 

FIG. 5. Field sweeping velocity dependence of the normalized gap. The noise is applied in 
each direction for (a) x-direction, (b) y-direction, and (c) z-direction. The solid line is the slow 
fluctuation limit of the noise. The noise has all components for (d). 

FIG. 6. Field sweeping velocity dependence of the normalized gap for the slow fluctuation limit 
of the noise. 
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Fig. 1 M. Nishino, et al. 




Fig. 3 M. Nishino, et al. 
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